Effects of the quark field on the ghost propagator of Lattice Landau Gauge QCD 
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Infrared features of the ghost propagator of color diagonal and color antisymmetric ghost propaga- 
tor of quenched SU(2) and quenched SU(3) are compared with those of unquenched Kogut-Susskind 
fermion SU(3) lattice Landau gauge. 

We compare 1) the fluctuation of the ghost propagator, 2) the ghost condensate parameter v of 
the local composite operator (LCO) approach and 3) the Binder cumulant of color anti-symmetric 
ghost propagator between quenched and unquenched configurations. 

The color diagonal SU(3) ghost dressing function of unquenched configurations has weaker singu- 
larity than the quenched configurations. In both cases fluctuations become large in g < 0.5GeV. 

The ghost condensate parameter v in the ghost propagator of the unquenched MILCc configuration 
samples is 0.002 ~ 0.04GeV^ while that of the SU(2) PT samples is consistent with 0. 

The Binder cumulant defined as U{q) = 1 — where (j){q) is the color anti-symmetric ghost 

propagator measured by the sample average of gauge flxed conflgurations via parallel tempering 
method becomes ~ 4/9 in all the momentum region. The Binder cumulant of the color antisymmetric 
ghost propagator of quenched SU(2) can be explained by the 3-d gaussian distribution, but that of 
the unquenched MILCc deviates slightly from that of the 8-dimensional gaussian distribution. 

The stronger singularity and large fluctuation in the quenched configuration could be the cause 
of the deviation of the Kugo-Ojima confinement parameter c from 1, and the presence of ordering 
in the ghost propagator of unquenched configurations makes it closer to 1. 

PACS numbers: 12.38.Gc, 12. 38. Aw, ll.15.Ha, ll.15.Tk 



I. INTRODUCTION 

Infrared features of the ghost propagator are important 
in the analysis of color confinement mechanism and the 
running coupling. Kugo and Ojima[3| considered the two 
point function connected by the ghost propagator and 
expressed the confinement criterion as 



i(0) = 1 - c = — = ^ 



(1) 



at the renormalization point ^ = 0f3|. Here Zi and Zi 
are the vertex renormalization factor of the triple gluon 
vertex and the ghost anti-ghost gluon vertex, respectively 
and Z3 and Z^ are the wave function renormalization 
factor of the gluon and the ghost, respectively. 

If Zi is finite, divergence of Z^ is a sufficient condition 
of the color confinement. The lattice data suggest that 
Z3 is infrared divergent, but its singularity is not strong 
enough to hinder the running coupling measured as 



gl Z{e)G{q^f 



in 



as{Kuv)q 



-2{an+2aG) 



(2) 
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approach zero in the infrared^. Here Z{q^) and G{q^) 
are the gluon dressing function and the ghost dressing 
function, respectively. The same observation is reported 

The ghost propagator in the infrared region was inves- 
tigated by several authors. Common findings are that 
it is more singular than and that in the infrared re- 
gion its statistical fluctuation is large probably due to 
presence of Gribov copies0,0,[3,BI3- quenched 
4g4 g^^^ gg4 gu(|3^ lattice simulation, the color di- 
agonal ghost propagator showed singularity of g^^-ac 
with ac ~ 0.25. In the Dyson- Schwinger (DS) ap- 
proach, the infrared power behavior of the ghost prop- 
agator and the gluon propagator j^ave the rela- 
tion 2aG -\- a£) — and the lattice data are consistent 
with this ansatz in q > IGeV region. As the magni- 
tude of an, Dyson-Schwinger approach^lQj and Langevin 
approachll jl p redicts -0.59, while the lattice data and DS 
approaches |l2l [Tsj predicts -0.5. If is smaller than - 
0.5 the gluon propagator in the infrared vanishes and the 
Gribov-Zwanziger's conjecture on the color confinement 
of the gluon becomes satisfied. Recent detailed analysis 
of the finite size effect in the lattice confirms that infrared 
limit of —ao in the DS app roach k — 0.5 is compatible 
with the lattice data[T^ll5l|. The relation 2aG + an ^ 
suggests presence of an infrared fixed point 10] . The in- 
frared finite quark wavefunction renormalization Z^ of 
unquenched simulation^^ also suggests that the running 
coupling is not infrared vanishing. 
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We cannot measure the ghost propagator at zero mo- 
mentum, since we evaluate it with the condition that it 
is zero-mode-less. Thus the infrared power fitted at finite 
lattice momentum ac cannot predict the power behavior 
of the ghost propagator near momentum i.e. the index 

K. 

In Q, we observed that the Kugo-Ojima confinement 
criterion is satisfied in the unquenched simulation but 
not in the quenched simulation of lattice sizes up to 
56** . In order to study the role of fermion in the color 
confinement, we consider the BRST(Becchi-Rouet-Stora- 
Tyutin) quartet mechanismjp. [l7| . 

In the BRST formulation[l| , unphysical degrees of free- 
dom are confined by the quartet mechanism. In the pure 
QCD in the Landau gauge, one can construct BRST 
quartet as 

-> D^{A)c ^ 

Af,c -> Df,{A)cc-A^B 0. 

Here the arrow implies the BRST transformation S b and 
B is the Nakanishi-Lautrup auxiliary field. The trans- 
verse gluon state is a BRST parent state of a daughter 
state Dfj_{A)c and the state with opposite ghost number 
of the D^{A)c i.e. Ac becomes a parent state, whose 
daughter and the above three states construct a quartet. 

Inclusion of the fermion field ip allows to construct an- 
other BRST quartet as 

ip —ipc 

Tpc —^cc — tpB — > 0. 

The Dirac fermion state V' is a BRST parent state of V'c 
and the state with opposite ghost number state of ipc is 
■00, which becomes a parent state of the BRST partner 
that construct a quartet. 

Inclusion of fermion gives more restiction on the de- 
grees of freedom of the ghost and it may change the fluc- 
tuation of the ghost propagator. 

Another current problem concerning the ghost propa- 
gator is the possibility of the ghost condensates. In the 
lattice Landau gauge QCD simulation, presence of 
condensates was suggested 0, IE IS Since A-^ is not 

BRST invariant, a mixed condensate i.e. a combination 
with ghost condensates 

/ {tra/HilA^A^' ~ ^tCC])d^x (3) 

was proposed|3ll3 ^ on-shell BRST invariant, i.e. in- 
variant for the B field that satisfies 

= --d^.A'"' + i^f'^'^c^c". (4) 

Here G/H is the subset of gauge fixed configuration, and 
^ is the gauge fixing parameter. The Landau gauge ^ = 
is regarded as a specific limit of the Curci-Ferrari gauge. 



In recent studies, the space-time average of the vacuum 
expectation value 

is claimed to have gauge invariant meaning 1 1 9l I2H . In 
the Landau gauge QCD, the Faddeev-Popov(FP) gauge- 
fixing action is 

Sfp = B-d^A';^ + ird^D^'^'c', (5) 

where the last term d ^"'^ C'' , where Df = (5"'' 9^ + 
gr'AI,. 

In analytical studies in the Curci-Ferrari gauge, pres- 
ence of the ghost condensate (/"'"^c''c'^) was discussed 
as the Overhauser effect in contrast to the {/'^'"^c^c'^) or 
(^fabc^b^c^ which are regarded as the BCS effect^. 

Since the Landau gauge is a specific limit of the Curci- 
Ferrari gauge, it is of interest to study the ghost prop- 
agator. In [g, we observed that in the SU(2) (3 = 2.1 
16'* lattice, the expectation value of color off-diagonal 
ghost propagator (e°'"^c''c'^) is consistent with but the 
standard deviation of the color-diagonal ghost propaga- 
tor has the momentum dependence of a{G°'°'{q)) oc q~'^. 
The investigation was extended by |2^ and this fluctu- 
ation was confirmed and although the expectation value 
of 4'"'{q) — e'^'^'^c^cf is consistent with 0, the expectation 
value of its absolute value was shown to behave 

as and not zero. We extend this approach to un- 
quenched MILC configurations. 

In [23, the ghost condensate parameter v and the 
Binder cumulant|24j of the color anti-symmetric ghost 
propagator was measured. In the Binder cumulant of an 
order parameter, renormalization factors cancel and one 
can extract the fixed point in the continuum limit by a 
suitable extrapolation. 

In the Zwanziger's Lagrangianj25j, the color anti- 
symmetric ghost field 4>''i^{x) leads to the mass gap equa- 
tion 

r'^{A'^^{x)cp'^{x)) = ^^^^^^^^^ (6) 

where 7^ is the mass dimension two Gribov mass 
parameter p^. It is not evident that the Zwanziger's La- 
grangian expresses the effective theory of the lattice Lan- 
dau gauge QCD, but a naly tical calculation of the ghost 
propagator in two loop|2q. and the local composite op- 
erator approach|22 |23| suggest hints for solving entan- 
glements in the confinement problem. 

In this paper we study the ghost propagator of 
quenched SU(2) (3 = 2.2 lattice gauge fixed to the 
Landau gauge via parallel tempering(PT) methodp^l 
and investigate the Binder cumulant. We extend 
the study to unquenched SU(3), using the MILCc 
configurations ^2^ . 

Organization of the paper is as follows. In sect. II, 
we show definitions of the color diagonal and color anti- 
symmetric ghost propagator on the lattice. In sect. HI, 



3 



fluctuation of the ghost propagator of quenched and un- 
quenched configurations are compared. In sect. IV, we 
compare the parameter v of the ghost condensates from 
the color anti-symmetric ghost propagators of quenched 
SU(2) PT configurations and unquenched MILCc config- 
urations. In sect.V, the Binder cumulant of the color 
anti-symmetric ghost propagator of the quenched SU(2) 
and unquenched SU(3) are compared. Summary and dis- 
cussion are given in sect. VI. 



II. THE GHOST PROPAGATOR 

The ghost propagator Da'^^{q^) and the ghost dressing 
fmiction G'^''{q^) is defined by the Fourier transform(FT) 
of the expectation value of the inverse Faddeev-Popov 
operator M — —d^V^ 

FT[DG^\x,y)] = FT{tr{A'^^{iM[U])-'y,yA')), 

= DG^\q^) = (7) 

where antihermitian SU(3) generator A'' is normalized as 
trA^U'' = 5''^. 
We measure 



DG^^q^) ^ {tr{k^q\M[U]-^\K\)) 



using the source vector |A"g) 



1 



A^e 



(8) 



We select 



the momentum to be directed along the diagonal of 
the lattice momentum space. 

In the approach of calculating the fourier transform of 
Sq{x))^, compensation of hypercubic artefacts was 
necessary, but in our method the artifact-free momenta 
are selected and the translation invariance is fully utilized 
to improve the statistics. 

The Faddeev-Popov operator M[U] = ~dD[U] is de- 
fined with use of the covariant derivative as 



(9) 



where d^cj) = (j){x + fi) — 0(x), (j) = + A*) + 0(a;)). In 

the tZ-linear version, {A^^^ = \{Ux^fj. - C^J,^)|tWp where 
\trip means the traceless part) S{Ux.fj.)Bx,f_i is defined as 



S{Ua:,fj.)Bx^^ - -{- 



and in the log— [/ version, (C/x,, 



Ax.fj 



2tanh(Ax,;./2) 



Bx,fj.j 



(10) 



(11) 



where Ax,f_, = adjAx,fj.^. 

In , the Faddeev-Popov operator was parametrized 



as 



Mifix, y) = S'^S{x, y) - f'^'^A^ix, y) (12) 



The authors decomposed the inverse matrix Dg^'^{x, y) = 
{M-')''%x,y) into DGl'{x,y) and Dgo'(^,?/) i.e. the 
component containing even number of A odd number of 
A, respectively. They derived the ghost propagator from 

gbc 

— DGf{x,y), and the color antisymmetric ghost 

1 

propagator by multiplying S^^A to the color antisym- 
metric ghost propagator DG^^{x,y) which contains per- 
turbation series of even numbers of A. 

We do not adopt this procedure, but derive directly 
the color anti-symmetric ghost propagators by the con- 
jugate gradient method. The convergence condition on 
the series is set to less than a few % in the I2 norm. 

We define M — —d^D^^ and solve 



A sinq • x 



and 



-a^Z?^/,^(x) = ^A^osq.x. 



(13) 



(14) 



Then we calculate the overlap to get the color diagonal 
ghost propagator 



DG{q) = 



1 1 

N}~IV 



(J"''((A° cosq-xl/^^x)) + (A°sinq-x|/,^(x))) 



(15) 



and color anti-symmetric ghost propagator 

xr'^iiA^ cosq • x|/,^(x)) - (A" sinq • x|/,^(x))) 



(16) 



where TV = 2 for SU(2) and 6 for SU(3). 



HI. FLUCTUATION OF THE GHOST 
PROPAGATOR 

We present in the following subsections square and the 
absolute value of the color antisymmetric ghost propa- 
gators of the quenched SU(2) /3 = 2.2 16^ lattice, and 
compare the corresponding values of SU(2) Landau gauge 
QCD of larger samples 23]. We measure also ghost prop- 
agators of quenched SU(3) (3 = 6.45 56^ lattice and 
those of unquenched MILCc with lattice size 20^ x 64 
and MILC/ with lattice size 28"^ x 96. We present square 
and the absolute value of color antisymmetric ghost prop- 
agator of MILCc. 



A. Quenched SU(2) 

We select momenta q following the cylinder cut, and in 
the case of unquenched SU(3) 20'^ x 64 lattice calculation. 
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it takes about 260 iterations in the q = 0.2GeV region but 
several iterations in the q = 4GeV region. The average of 
color anti-symmetric ghost propagator (j)'^{q) is consistent 
with but the average of its square (t>'^{qY has a non- 
vanishing value. We define 



4>{qf 



(17) 
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FIG. 1; Log of color anti-symmetric ghost propagator squared 
logi(,[<^(g)^] as the function of q{GeV). f3 = 2.2, 16" PT gauge 
fixing. (Color online) 

The log of ${q)'^ of /? = 2.2 16** lattice gauge fixed by 
the PT method (67 samples) is shown in Fig^ The cor- 
responding log-log plots (Fig[21) is to be compared with 
that of the ghost propagator Dc{q) (Fig|2J). 
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FIG. 2: Log of the color anti-symmetric ghost propagator 
squared logj^Q[(/f)(g)'^] as the function of logjQ[g(GeV)]. /3 — 2.2, 
16'' PT gauge fixing. (Color online) 



The infrared singularity of the standard deviation of 
the color anti-symmetric ghost dressing function and 
the color diagonal ghost dressing function are and 
q~'^'^, respectively. 



B. Quenched SU(3) 

In FIG 0] we show the color diagonal ghost propaga- 
tor of quenched SU(3) with (3 = 6.4 (1/a = 3.66GeV) 
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FIG. 3: Log of the ghost propagator \og-^Q[DG{q)] as the func- 
tion of logio[g(GeV)]. 13 = 2.2, 16'' PT gauge fixing. (Color 
online) 
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FIG. 4: The ghost propagator as the function of the mo- 
mentum g(GeV). /3 = 6.45, 56'(stars) and (3 = 6.4, 56'(filled 
diamonds) in the log U definition. Solid line is the pQCD fit 

in MOM schemeQ. (Color online) 



and (3 = 6.45 (1/a = 3.8697GeV) on 56-* lattice. The 
corresponding ghost dressing function are in FIG[S1 

The standard deviation of the color-diagonal ghost 
propagator oi (3 = 6.45 multiplied by {qaY is almost 
constant in g > IGeV region but in the q < O.SGeV re- 
gion it is enhanced as compared to the value at q > IGcV 
region. 
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FIG. 5: The ghost dressing function as the function of the mo- 
mentum g(GeV). P = 6.45, 56'(stars) and /? = 6.4, 56'(filled 
diamonds) in the log —U definition. (Color online) 
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The log — log plot of the standard deviation of the 
color-diagonal ghost dressing function of /3 — 6.45 in the 
q < IGeV region behaves as 



a(G(g)) oc (7-2.8(1)^ 



C. Unquenched SU(3) 



(18) 



In we showed lattice results of color diagonal ghost 
dressing function of unquenched JLQCD/CP-PACS and 
MILC. In these simulations the length of the time axis is 
longer than the spacial axes and the ghost propagator of 
low momentum region is extended. In FIGElthe log- log 
plot of the ghost dressing function of the MILC / (3imp = 
7.09 on 28^ x 96 lattice and that of quenched /3 = 6.45 
on 56"* lattice are shown. We observe suppression of the 
ghost pro pag ator in the infrared region in the asymmetric 
lattice^, lig- Systematic deviation of ghost propagator 
and gluon propagator of asymmetric lattice from those 
of symmetric lattice is recently confirmed in the large 
3-dimensional SU(2) lattice|30|]. The suppression in the 
infrared of the unquenched data may not be due to the 
presence of quarks but due to the geometry of the lattice. 






u 
o 



5 



-0.5 







0.5 



Or 

^ -0.5- 
-1- 

,r -1.5- 

S -2- 
^-2.5- 

o 

tji -3 - 
o 

^^ -3.5- 



-0.8 -0.6 -0.4 

Logio [q(GeV) ] 



-0.2 



FIG. 7; Log of the standard deviation logj^Q a{G(q)) as a func- 
tion of logiQ (jr(GeV) of MILC/ /3imp = 7.09 (upper points) and 
that of quenched /? — 6.45 56* (lower points). (Color online) 



The color anti-symmetric ghost propagator of MILCc 
(21 samples) is shown in Fig|S| By comparing with 
FIG^ we observe decrease of the slope. 
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FIG. 8: Log of the color anti-symmetric ghost propagator 
squared logjQ [</>(g)'^] as the function of g(GeV). Pimp = 6.83 
and 6.76, 20^ x 64 MILCc. (Color online) 
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FIG. 6: Log of the ghost dressing function log^Q G{q) as a 
function of logj^Q (j'(GeV) of MILC/ /3imp = 7.09 (diamonds) 
and that of quenched f3 — 6.45 56* (stars). (Color online) 

There are difference in the momentum dependence of 
the standard deviation of the color diagonal ghost dress- 
ing function of unquenched MILC/, (3imp — 7.09 on 
28'^ X 96 lattice and the quenched (3 = 6.45 on 56"* lat- 
tice in the infrared region as shown in FIG[71 Since the 
sample size is different, the absolute value of the stan- 
dard deviation is not meaningful, but the strength of the 
fluctuation defind by the slope influences the infrared be- 
havior of the running coupling etc. 

The momentum dependence of the standard deviation 
of the color diagonal ghost dressing function of MILC / is 
dramatically less singular than that of quenched config- 
uration. We observed 

a(G(g))(xg-i-iW. (19) 



IV. THE GHOST CONDENSATE 

In the case of SU(2), the ghost condensate appears in 
the color anti-symmetric ghost propagator DQ'^{q) re- 
lated to 0°((?) through 

£ahc £ahc 

= ~^^DG''{q') = -i^Dc'^iq) (20) 

and 

DG'':iq)^^'^^^e^'^, (21) 
q* + 

where e^^ is an anti-symmetric tensor, i.e. when a = 3, 6 
and c = 1,2. In general, we parametrize the average of 
I </-"(?) I as 

^ a 



6 



Here L is the lattice size and the parameter r/L^ is 
the correction from the finite size effect. 



A. Quenched SU(2) 

In (23I, the fitting parameter r of (color anti- 

symmetric ghost propagator) on the lattice was derived 
from 



7^^l'^°('?)l = ^' (23) 
3 ^ cos{nq/L) 



in which q = 0,1, ■ ■ ■ L. 

Our fit of — ^f* ^'^'^j — - of PT samples using r 
cos[TTqa/ L) 

10.13, z = 4.215 is shown in FIGEI 



^ 
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B. Unquenched SU(3) 

As in the SU(2) PT samples, we performed the fit of 
the parameter v for the MILCc samples. We first fit 
the log of |0°(g)|(color anti-symmetric ghost propaga- 
tor) devided by cos{nqa/ L) using L — v2(?^x64 and 
obtained r = 40.5, z = 3.75, as shown in FIGHH The 
parameter r and z for the fit of cos^(7r^/L) are 

r = 36.5, z — 7.5. 



to 



-0.75-0.5-0.25 0.25 0.5 0.75 

Logioq(GeV) 

FIG. 11: Log of |(/>(g)|/ cos{TTqa/ L) (color anti-symmetric 
ghost propagator) as the function of logjQ(g(GeV)) of MILCc 
samples. (Color online) 

Our fit of |(/'((7)| ignoring two lowest momentum points 
and using r = 40.51 gives v = 0.0020GeV^, which is small 
but positive. When the two lowest momentum points are 
included, v decreases to -0.0005 but x^/d-o.f. increases. 
The former fit is shown in FIGI12I 




FIG. 9: Log of 10" (g) I (color anti-symmetric ghost propaga- 
tor) devided by cos(7rga/L) as the function of logj^Q(q(GeV)) 
of SU(2) PT samples. (Color online) 

Using r = 10.13, L = 16, the fitting parameter v of 
|0(g)| is found to be -0.002GeV^ and is consistent with 
0. The fit is shown in FIGCUl 
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FIG. 10: Log of the absolute value |(;/>°(g)| (color anti- 
symmetric ghost propagator) as the function of logjQ(g(GeV)) 
of SU(2) PT samples. (Color online) 
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FIG. 12: Log of the absolute value 10(g) | (color anti- 
symmetric ghost propagator) as the function of logj^Q(q(GeV)) 
of MILCc samples. (Color online) 

We fitted also logjQ0^((?), where 

^ a 

The fit with r ^ 40.5, v = 0.035GeV2 is shown in 
FIGini The fit with r = 36.5, v = 0.041GeV2 is not 
distinguished from this figure. 
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FIG. 13: Log of the 4>{q)^ (color anti-symmetric ghost prop- 
agator squared) as the function of logjQ(g(GeV)) of MlLCc 
samples. (Color online) 



V. BINDER CUMULANT 

Two decades ago Binder 24] showed cumulants of the 
order parameter yields non-trivial fixed-point values. 
The theory was applied to the Ising model in which the 
magnetization M is the order parameter jsH Is^ and the 
cumulant was defined as 



B 



2^ (M2)2' 



(25) 



When the distribution of M is given by the 1 dimensional 
gaussian distribution, one finds 



(26) 



and B becomes 0. 

In SU(2) and SU(3) lattice QCD, deconfinement phase 
transition was studied by measuring 



(P2)2 



(27) 



using the Polyakov line data P as the order parameter p3a 

Since the color anti-symmetric ghost propagator could 
be an order parameter of the system, the authors of 
considered its Binder cumulant defined as 



U{q) - 1 - 



3(0(9)2)2 • 



(28) 



We measure 



—y\- 



x^— ((A''cosq-x|7W"i|A'=sinq-x) 
-(A''sinq-x|X-i|A=cosq-x))]2 



a 

{{K^ cos q • x| A^-i I A'^ sin q • x) 



/ 



V 



-(A''sinq-x|7W-i|A'=cosq-x))]2)2. (30) 

In arbitrary d-dimcnsional space, corresponding expec- 
tation value for d-dimensional gaussian distribution be- 
comes 



rf + 2 



(31) 



Thus a natural extension to d-dimensional vector variable 
is 



U{q) 



d + 2 



(32) 



which becomes in the system with gaussian distribution 
whose symmetry is not broken. 

When the symmetry of the system is broken, as in the 

Ising model at the temperature, the ratio of {(f)[q)'^)'^ 

~ 2 
and {(f){qY) becomes 1 and Uiq) becomes . It corre- 
ct 

sponds to the temperature fixed point. 



A. Quenched SU(2) 

We measure the Binder cumulant U{q) of the quenched 
SU(2) 16^ 13 = 2.2, a = 1.07GeV-i configurations (67 
samples) produced by the PT Landau gauge fixing and 
the corresponding first copy[^ The q dependence 

of the U{q) of PT gauge fixed samples and the first 
copies are shown in FIG. El The infrared fluctuation is 
large in the first copy but it is reduced in the PT gauge 
fixed samples. It implies a Gribov copy effect in the in- 
frared region The average over q > O.SGeV becomes 
U{q) = 0.45(2). This value is comparable to that of 
obtained by 10000 samples using symmetric momentum 
91 = 92 = 93 = (74 7^ 0. In [23, the value U{q) be- 
tween and 2/3 was interpreted as a system deviating 
from the gaussian distribution. However, since (f>{q) is 
a 3-dimensional vector, it would not be appropriate to 
treat it as an 1 dimensional object. The value 0.45 is 
very close to 



Uiq) ^ 1 - 



rf + 2 _ 4 
3d ~ 9' 



or 



(29) corresponding to the 3-d gaussian distribution. 



(33) 



(34) 
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FIG. 14: The momentum dependence of Binder cumulant 
U{q) of SU(2), /3 = 2.2, a = 1.07GeV"^ of PT samples and 
first copy samples, (color online) 



B. Unquenched SU(3) 

We measured the Binder cumulant of the color anti- 
symmetric ghost propagator of MILCc. We observed 
qualitatively different features from quenched SU(2). An 
average of 9 Pimp = 6.76 samples and 12 Pij^p = 6.83 
samples of MILCc is shown in Fig^| When the (f> is dis- 
tributed as a Gaussian vector in eight dimensional space, 
U{q) = 1 - = 0.58. Data of FigUS suggests that 
U{q) is slightly larger than the 0.58, and that the shift 
from Gaussian distribution of Pimp = 6.76 samples with 
light bare quark mass niQ = 11.5MeT^ is larger than that 
of Pimp = 6.83 samples with heavier bare quark mass 
mo = 65.7Mey. 

A qualitative difference of unquenched SU(3) fFig [TB|) 
from quenched SU(2) fFig ll4|l is the smallness of the fluc- 
tuation at the lowest and next to the lowest momentum 
point [q = (0,0,0, 1) and (0,0,0,2)]. Relatively large fluc- 
tuation exists when one of the spacial components of q is 
1 and other components are 0. The difference from the 
quenched SU(2) could be due to the improvement in the 
Asqtad action used in the MILCc and/or the presence of 
dynamical fermions. 



4j 
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FIG. 15: The momentum dependence of Binder cumulant 
U{q) of unquenched SU(3), l3imp = 6.83 and 6.76, a — 
1.64GeV-^ MILCc. (color online) 



VI. SUMMARY AND DISCUSSION 

We presented color diagonal ghost propagator of 
quenched P — 6.45 56^ lattice and those of unquenched 
MILCc 20^ X 64 and MILC/ 28^ x 96 configurations. The 
momentum dependence of standard deviation of the color 
diagonal ghost dressing function of the unquenched con- 
figurations is less singular than that of the quenched con- 
figurations. The standard deviation and the mean value 
of statistical distribution is important for determining 
the nature of the ensemble. 

The ghost pair creation operator in the BCS channel 
is expected to behave as the order parameter, and in the 
Landau gauge, where ghost pair creation is absent, the 
ghost anti-ghost pair creation in the Overhauser channel 
was speculated to become an order parameter. The pa- 
rameter V of LCO approach that characterize the ghost 
condensate was compatible with in the SU(2) PT sam- 
ples. In the unquenched SU(3) MILCc samples, we found 
small but positive value of v. Uncertainty on v comes 
mainly from that of r, where finite size effect is crucial. 
We need to extend the analysis to larger lattices to get 
the definite conclusion. 

We showed that the Binder cumulant which mea- 
sures the fluctuation of the ghost propagator differ be- 
tween quenched and unquenched conflgurations. We con- 
flrmed that the Binder cumulant U{q) of the color anti- 
symmetric ghost propagator of SU(2) obtained by 10000 
samples 23] U{q) ^ | is consistent with that obtained 
by using PT gauge flxed samples. In 3-dimensional sys- 
tem, this data can be interpreted as U{q) defined as 
ea. ()34|l ~ 0, i.e. the color symmetry is not broken and 
that the system is in the random phase. When the sys- 
tem is ordered, a certain direction in the color space will 
be selected and the Binder cumulant would deviate from 
the value expected by the gaussian distribution. The 
data of quenched SU(2) P — 2.2 do not show this ten- 
dency but the unquenched SU(3) show deviation from the 
gaussian distribution. Whether it implies the precourser 
of the ghost condensation in the unquenched QCD is not 
evident. It would be interesting to extend the analysis to 
finite temperature and study qualitative differences. 

The ghost condensates and the condensates are ex- 
pected to be related by the on-shell BRST symmetry. 
The observables of lattice Landau gauge in l-3GeV re- 
gion suggests presence of condensates. The larger 
standard deviation of the SU(3) quenched ghost propa- 
gator as compared to the unquenched ghost propagator 
may imply that the ghost propagator is more random 
in the quenched samples. The fluctuation of the ghost 
propagator could be the main cause of suppression of the 
running coupling in the infrared and saturation of the 
Kugo-Ojima parameter c at about 0.8 in the quenched 
approximation. It is likely that the fermion field reduces 
the fluctuation of the color diagonal ghost propagator in 
the infrared, and renders the Kugo-Ojima parameter c 
close to 1. 

We think infrared suppression of the running coupling 
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of unquenched SU(3) measured by eq.(|51) presented in Q 
is a finite size effect. In the process of measuring the 
ghost propagator for the running couphng, we observed 
large fluctuations of the norm and random orientation 
of the vector in adjoint color space i.e. weakening of 
the color diagonal structure of the ghost propagator in 
the infrared. Concerning the fixed points of the running 
coupling, Wilsonls^ noted in 1971 that the renormaliza- 
tion group flow of the coupling could approach limit cy- 
cles which are more elaborate than simple isolated fixed 
points. A possibility of complicated fractal structure in 
fixed points was discussed in The Zamolodchikov's 
c-theorem in two dimensional conformal field theorv,:37.J, 
however, excludes the limit cycle structure of the infrared 
fixed points. In four dimensional QCD, the situation is 



obscure'Sq. To clarify the nature of the infrared fixed 
points, it is necessary to investigate the continuum limit 
of the lattice Landau gauge QCD via systematic studies 
of finite size effects and the Gribov copy effects. 
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